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We experimentally demonstrate optimal entanglement distillation from two forms of two-qubit 
mixed states under local filtering operations according to the constructive method intruduced by 
F. Verstraete et al. [Phys. Rev. A 64, 010101(R) (2001)]. In principle, our set-up can be easily 
applied to distilling entanglement from arbitrary two-qubit partially mixed states. We also test the 
violation of the Clauser-Horne-Shinmony-Holt (CHSH) inequality for the distilled state from the 
first form of mixed state to show its " hidden non- locality" . 
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Entanglement plays a key role in quantum information processing, such as quantum teleportation 0, efficient 
quantum computation and entangled-assisted quantum cryptography Q. In general, these applications require 
maximally entangled quantum states. However, owing to decoherence and dissipation, practical states are normally 
less entangled or partially mixed. To c ®p e with this problem, entanglement concentration is essential Q and various 
schemes have been proposed 0, S 0, U Experimentally Kwiat et al. implemented entanglement distillation from 
non-maximally entangled pure states and one special kind of two-qubit mixed states using partial polarizers ||| ■ After 
that, ent ang lement concentration from maximally entangled mixed states (MEMS) has also been realized using similar 
i method jlOj . 

' Generally speaking, there are two types of concentration protocols: those involving collective operations performed 

■ on many copies of the state and those working on individual copies. The latter case is of special interest in practical 

■ protocols, because usually there might be technologic difficulties to implement collective operations. For this case, 
although some experiments on certain kind of two-qubit mixed states have been experimentally demonstrated j^, llC| , 
it is still interesting to explore what is the best we can do to distill entanglement from arbitrary two-qubit mixed 
states, only allowing local operations and classical communication (LOCC) to be performed on each copy separately. 
In Ref. p], Kent et al. proved that the best state one can obtain from general two-qubit mixed states is a Bell 
diagonal state. After that, F. Verstraete et al. constructively gave out the optimal local filtering operations for 
distilling entanglement from an arbitrary two-qubit mixed state In this paper, we experimentally demonstrate 
optimal entanglement distillation from two forms of two-qubit mixed states which frequently exist in the real world. 
The optimal local filtering operations are explicitly calculated according to the method provided in Ref. p~l| . And in 
principle, our set-up can be applied to entanglement distillation of arbitrary two-qubit mixed states. 

For general two-qubit mixed states, consider two separated parts, Alice and Bob, who each control one subsystem 
H ' and are only allowed to carry out local operations and classical communications. Specifically, Alice and Bob are only 
permitted to perform local unitary transformations and local filterings. Then what we are concerned about is whether 
•rH . they can increase the entanglement of the system and to what extent they can do that under LOCC. This problem 
1 is very interesting because any real world quantum communication channel is imperfect and generally, Alice and Bob 
are not able to share perfect maximally entangled states directly. 

Theoretically, by introducing the real and linear parametrization, a two-qubit density matrix can be represented 
as P = lEi,3=0'^»J (7 > ® a j> wnere cr is the 2 x 2 identity and o-i,a 2 ,o- 3 are usual Pauli matrices E|. LOCC of 
the type p ~ {A <£> B)p(A ® B)^ corresponds to left and right multiplication by a proper orthochronous Lorentz 
transformation (POLT) , followed by normalization in this i?-picture. The concentration protocol, that is optimal in 
the sense that it produces a state of maximum entanglement of formation (EoF) Q , has been obtained in [Tll | for two 
cases. If R = [Rij] (Rij — Tr(p(ai ® 0j))) is diagonalizable by POLT, a Bell diagonal mixed state can be extracted 
from the input state with a finite probability, which has the maximum EoF and the maximum possible violation 
of the Clauser-Horne-Shimony-Holt (CHSH) version of inequality [T3, Il3[. When R is not diagonalizable by POLT, 
the initial state can be quasi-distillable. In the extreme case, the input state can be asymptotically transformed 
into a Bell diagonal state with lower rank, while the success probability becomes infinitesimally close to zero. This 
method can also be used for entanglement concentration from MEMS since the method in |l(J is a particular case of 
it. 
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Experimentally, to demonstrate the optimal entanglement distillation protocol, we concentrate on the mixed states 
that can be diagonalized by POLT. We specially devise two forms of two-qubit mixed states which can be easily 
prepared and distilled under fewer local operations. The experiment set-up to investigate entanglement distillation is 
shown in Fig. 1. A 0.59 mm thick /3-barium borate (BBO) crystal arranged in the Kwiat type configuration is 
pumped by a 351.1 nm laser beam produced by an Ar + laser. In the spontaneous parametric down-conversion (SPDC) 
process, a nonmaximally entangled state a\HH) + b\VV) ( H and V represent horizontal and vertical polarization of 
the photons respectively) is produced. In the experiment, for simplicity, we adjust the half wave plate (HWP) and 
tiltable quarter wave plate (QWP) in the pump path to choose the real numbers a and b [T^j. The experiment can 
also be applied to the situation where a and b are complex numbers. After the BBO crystal, the photon pairs pass 
through two phase-damping channels in a particular basis. In our scheme, each phase-damping channel consists of 
one quartz plate with its optical axis rotated by a certain angular to the vertical axis. For the first form of mixed 
state, we select two same phase-damping channels in {if + V, H — V} basis with the corresponding super-operators 
{^/l — pi, s /p<j x }a ® j-y/l — pi, y/p<7 x }B, where p is connected with the thickness of the quartz and the bandwidth of 
the interference filter [lj|. Then the mixed state has the form 

/ (p- 1) 2 + b 2 {2p- 1) ab((p- l) 2 +p 2 ) 

p-p 2 2abp(l-p) 

2abp(l-p) p-p 2 

\ ab((p- l) 2 +p 2 ) p 2 -b 2 {2p-l) 



Pi 



(1) 



In the experiment, we choose a — 0.23, b — 0.97, and p = 0.013 (a ~ 1.5-mm-thick quartz plate). 

According to Ref. ^lj > the local operations for Alice and Bob are ( q 49 ) an< ^ ( 1 ) ( 49 ) respectively. 

For Alice, the filtering operation can be realized by inserting into her path microscope slides tilted about the vertical 
axis 73°. This configuration has a measured transmission of 92% for H polarization and 22% for V polarization. As 
for Bob, there are one local filtering and one single-qubit unitary operation. He just needs to place a half wave plate 
(HWP) after the microscope slides with the angle between the optical axis and the vertical axis set to 45°. We can 
also realize the local filtering operation using a Mach-Zehnder interferometer [l6| . but it requires too much about its 
stability and visibility. 

To obtain the density matrices of the distilled states, we use the technology of maximum likelihood tomography 0, 
\vi\ . By means of a quarter wave plate (QWP), HWP and polarizing beam splitter (PBS) in each arm, the polarization 
of each photon can be analyzed in an arbitrary basis. Then by combinations of 16 single-photon projections \H), 
\V), \H) + \ V), | H) + i\V) on each of the two-photon, we may derive the density matrices describing the states of 
the photon pairs (all the experimental density matrices in this paper can be found in 01)- Fig. 2 shows the density 
matrices before and after entanglement distillation. 

Here we use concurrence to characterize the entanglement between the two subsystems since the EoF of a two- 
qubit system is a convex monotonically increasing function of the concurrence hj. The concurrence before and 
after distillation are 0.248 ± 0.021 and 0.672 ± 0.044 respectively. The distilled state has the fidelity of 20] 82% 
with the theoretical Bell diagonal state. From the data we can see that by means of the local filtering operations we 
can effectively increase the entanglement between the two subsystems. The errors mainly stem from the imperfect 
preparation of the initial mixed state (the fidelity is 94%) and in the case of strong filtering, small changes in the 
initial state may have a large impact on the final state. 

We also show the "hidden non-locality" of the distilled state by its violation CHSH version of inequality. In the 
CHSH inequality, the proposed value \S\, a combination of four polarization correlation probabilities, should not 
be more than 2 for local hidden variables theory. If \S\ > 2, we can only use quantum mechanics to explain the 
correlations. Experimentally, the set-up for tomographic reconstruction can also be used for the measurement of the 
CHSH inequality. The analysis settings of QWPs and HWPs are determined by the tomographically measured density 
matrix of each state using the method introduced in Ref. [2l| . Following the normalization procedure j2^| , we obtain 
the value Sbefore = 1-853 ± 0.011 and S a fter = 2.175 ± 0.024, corresponding to the values of S before and after the 
entanglement distillation. The ideal value of S a ft er is 2.192. The experiment result accords well with the theory. 
Thus we obtain the violation of about 7 ex. In the process of entanglement distillation, Alice and Bob only employ 
LOCC operations. As pointed by Kwiat [9j, the non-locality demonstrations rely on "conditional probabilities": we 
use LOCC operations to select a subensemble with a small probability which can demonstrate non-local correlations. 
This idea is similar to entanglement distillation: we can just increase the entanglement of subensemble and can not 
do that for the whole ensemble [23| . 

Next consider distilling the entanglement from the other form of two-qubit mixed state. This time we let the 
nonmaximally entangled states pass through two same phase-damping channels in {H, V} basis. The correspond- 
ing super-operators are — pi, ^fpa z } a <8 — pi, \Jp~Vz\b- Each channel is a ~ 3 mm quartz which has a 
polarization-dependent optical path length difference about 40A (A = 702.2 nm). After the decoherence, the state has 
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the form 



Pn 



/ a 2 ab(-l + 2 P y 





\ab{-l + 2p) 2 b 2 
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For this kind of mixed states, the only local operation is just a unilateral local filtering f q ^ J f° r Bob. We consider 

two cases (a = 0.44, 0.52 and p — 0.063). The experiment results are shown in Fig. 3. Under this local filtering 
operation, the entanglement of the final states can be increased. For a = 0.44, the concurrence increases from 
0.552 ± 0.017 to 0.641 ± 0.022; the fidelities with the corresponding theoretical states are 98% and 96% respectively. 
For a = 0.52, the concurrence increases from 0.569 ± 0.017 to 0.666 ± 0.021; the fidelities with the corresponding 
theoretical states are 97% and 97% respectively. 

In conclusion, we have experimentally demonstrated optimal local filtering operations for two general forms of 
two-qubit mixed states that can be diagonalized by POLT. For the first form of mixed state, the local operations 
involve bilateral filtering; while for the second one, only unilateral filtering is needed. In fact, we can generalize this 
method to arbitrary two-qubit partially mixed states, because any nontrivial LOCC operations have the form 0, Il6| 

jUa y q j U' A ® Ub ^ g ^ ^ U' B , where Ua(b) aim ^A(B) denote local unitary operations for Alice (Bob) and 

7 is a scale factor in the range < 7 < 1 and < (xa(b) < 1- The main experiment difficulty lies in the preparation 
of arbitrary two-qubit mixed states. 

On the other hand, the imperfect preparation of initial mixed states is also a main reason for the deviations of 
our results from theoretical calculation. However, in practical schemes, one can directly choose the states derived 
from the tomographic measurment of the input as the initial mixed states and implement corresponding local filtering 
operations on them. Thus the experiment results will approach the theory with better accuracy. After the local 
filtering operations, Bell diagonal states can be obtained from the input states which have the maximum concurrence 
and the maximum possible violation of the CHSH inequality. We experimentally obtain the violation of the CHSH 
inequality about 7a for the distilled state from the first form of mixed state, which verifies "hidden non-locality" for 
this form of mixed state in the process of entanglement distillation. Recently it has been found that Bell diagonal 
states can also be used to implement the generalized tomographic quantum key distribution protocol |24j . The method 
used in this experiment is very simple and universal since it just needs LOCC operations and it can effectively increase 
the entanglement of arbitrary two-qubit partially mixed states. We believe it will be helpful in the exploration of 
various quantum information processing. 
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FIG. 1 Experimental set-up to implement entanglement distillation. The half wave plate, quarter wave plate and BBO 
crystal in the pump beam are used to prepare nonmaximally entangled states. In the mixed state preparation process, we 
use quartz plates as decoherence channels which can be rotated to introduce decoherence in a particular basis. The local 
operations are realized using microscope slides and HWPs. To distill entanglement from the first form of mixed state, the 
local operations involve bilateral filterings. For the second form of mixed state, only unilateral filtering is needed. The final 
QWP and HWP, together with polarizing beam splitter (PBS) in each arm, enable analysis of the poarization correlations 
in any basis, which can be applied to tomographic measurement and Bell inequality analysis. To detect the photon pairs, 
we use interference filters (bandwidth 4.62 nm), single-photon detectors and two-photon coincidence. 

FIG. 2 Experimental data for the first form of mixed state. Density matrices for the states before and after the entanglement 
distillation are shown in Figure (1) and (2). The fidelities with the corresponding theoretical states are 94% and 82%. 
FIG. 3 Experimental data for the second form of mixed state. Density matrices in the case a = 0.44(0.52) before and after 
the entanglement distillation are shown in Figure (1)((2)) and (l')((2')). The fidelities with the corresponding theoretical 
states are 98%(97%) and 96%(97%). 
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(1) Before distillation 



(2) After distillation 
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(1) Before distillation 



(1' ) After distillation 
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